Coherent back-scattering near the two-dimensional metal-insulator transition 
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We have studied corrections to conductivity due to the coherent backscattering in low-disordered 
two-dimensional electron systems in silicon for a range of electron densities including the vicinity 
of the metal-insulator transition, where the dramatic increase of the spin susceptibility has been 
observed earlier. We show that the corrections, which exist deeper in the metallic phase, weaken 
upon approaching to the transition and practically vanish at the critical density, thus suggesting 
that the localization is suppressed near and at the transition even in zero field. 

PACS numbers: 71.30.+h,73.20.Fz,73.20. Jc 



•4-* 



C 

O 

o : 

(N 
> 

oo 

o 

O 
O 

(N 
O 
-i— > 



T3 
C 

o 
o 



X 



The metallic conductivity and the metal-insulator 
transition (MIT) in two-dimensional (2D) systems have 
been a subject of intense studies for almost a decade 
Q. It is now well established that strong interactions 
between carriers are responsible for the metallic behav- 
ior, with quantitative agreement between the experi- 
ment and interaction-based theory 0, reported for 
a variety of dilute 2D systems. However, no quantita- 
tive theory is currently available in the most interest- 
ing regime: in the immediate vicinity of the transition, 
where electron-electron interactions cause strong changes 
in the system parameters including dramatic rise of the 
spin susceptibility 0, 0, IE possibly leading to a spon- 
taneous spin polarization at a finite electron density n x 
0, Q . One of the open questions is the fate of the co- 
herent backscattering (responsible for the weak localiza- 
tion) in the presence of strong interactions. In some 
experiments on p-GaAs/AlGaAs heterostructures, neg- 
ative magnetoresistance consistent with suppression of 
the coherent backscattering by a perpendicular magnetic 
field, B±, has been reported at carrier densities as low 
as 4 x 10 10 cm" 2 corresponding to r s as 10, which has 
been interpreted as evidence of the insulating character 
of the ground state (here r s = E c /Ep is the ratio 
between Coulomb and Fermi energies characterizing the 
strength of the interactions). However, in ultra-clean p- 
GaAs/AlGaAs heterostructures with much lower carrier 
densities (down to 7 x 10 9 cm" 2 ), the negative mag- 
netoresistance has been found to be negligible near the 
metal-insulator transition ^(J — less than 3% of its "non- 
interacting" value — suggesting that the effects of the 
coherent backscattering are suppressed in this strongly- 
interacting system. 

In this Letter, we report the first experimental studies 
of the weak-field magnetoconductance in high-mobility 
silicon metal-oxide-semiconductor field-effect transistors 
(MOSFETs) at temperatures down to about 40 mK and 
at electron densities down to the critical density for the 
MIT, n c = 0.82 x 10 11 cm" 2 . In this regime, spin sus- 
ceptibility sharply increases (and, possibly, even diverges 
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UlElH) indicating that prop- 



erties of the system are dramatically altered by the in- 
teractions. Previous studies of the coherent backscatter- 
ing in this system were restricted to electron densities 
n s > 1.5 x 10 11 cm" 2 11], which are outside of this 
interesting regime. We have found that relatively far 
from the transition, the corrections to the conductivity 
due to the coherent backscattering are of their "normal" 
(non-interacting) strength, in agreement with earlier re- 
ports 11] . However, once the transition is approached 
(8 = n s /n c — 1 < 0.1), the corrections become progres- 
sively weaker, and in the immediate vicinity of the critical 
density, they practically vanish. This suggests that the 
weak localization is already suppressed near and at the 
transition even in the absence of the perpendicular mag- 
netic field, in agreement with recent data on ultra-high 
mobility p-GaAs/AlGaAs heterostructures The ab- 
sence of localization at and just above the critical density 
may account for the existence of a flat separatrix between 
metallic and insulating states observed at n s = n c in 
Refs-HHH. 

We used three high-mobility samples from dif- 
ferent wafers which all displayed similar behavior. 
The data shown below were obtained on two sam- 
ples (NH50.3064.9N and NH50.4064.7N) from different 
wafers; peak electron mobilities in both samples were 
close to 3.3 x 10 4 cm 2 /Vs (at T = 0.1 K). The resis- 
tance was measured by a standard four-terminal tech- 
nique at very low frequencies (0.2 to 0.5 Hz) to mini- 
mize out-of-phase signal. The voltage between poten- 
tial probes was amplified by a factor of 10 3 by a ultra 
low-noise battery-operated SR-560 differential amplifier 
with a symmetric input and then measured by a SR- 
830 lock-in detector. There are two main experimen- 
tal difficulties associated with low-temperature transport 
measurements: the first is to minimize the contact re- 
sistance, which tends to grow very high at mK tem- 
peratures, and the second is to ensure that the carri- 
ers are not overheated relative to the bath temperature. 
The first problem was solved |14| by introducing gaps in 
the gate metalization which allowed for maintaining high 
electron density near the contacts even when the den- 
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sity in the main part of the sample was low. The gaps 
were narrow enough (< 100 nm) for the given gate-oxide 
thickness to provide a smoothly descending electrostatic 
potential from the high-density part to the low-density 
part. To minimize the second problem, we filtered out the 
external noise by using two sets of ir filters. The elec- 
tron gas was cooled by 16 copper wires approximately 
0.5 mm in diameter thermally coupled to the mixing 
chamber of a Kelvinox-100 dilution refrigerator. To ver- 
ify that the electrons were not overheated relative to the 
bath temperature, we studied the temperature depen- 
dence of the amplitude of the Shubnikov-de Haas oscilla- 
tions. In Fig. 1, we show p(B±) dependences recorded 
at "high" electron density, n s = 7.3 • 10 11 cm~ 2 (a), 
and at "low" density, n s = 1.17 • 10 11 cm" 2 (b). We 
have carefully checked that the data were taken in the 
linear regime of response. The temperature dependence 
of the normalized amplitude of the oscillations, AJAq, 
is shown in the insets together with the theoretical de- 
pendences calculated using Lifshitz-Kosevich formula 0] 
(here A — 4exp(—2ir 2 kBT£>/huj c ), uj c — eB±/mc is the 
cyclotron frequency, m is the effective mass, and To is the 
Dingle temperature). As seen from the figures, the theo- 
retical dependences coincide with the experimental data 
at both high and low electron densities down to the low- 
est accessed temperature rs 40 mK. This confirms that 
the electrons were indeed "cold" in our experiments. The 
power generated in the samples was maintained at about 
10 -15 W at low electron densities (n s < 1.3 ■ 10 11 cm~ 2 ) 
and at about 10 -14 W at higher electron densities. The 
measuring current ranged from 40 pA at the lowest den- 
sities to 2 nA at the highest ones. 

In Fig. 2 (a), we show dependences of the longitudi- 
nal magnetoconductance, a, on a perpendicular magnetic 
field recorded at T — 42 mK for the range of electron 
densities on the metallic side of the metal-insulator tran- 
sition. All curves are symmetric about B± = 0. At 
higher n s (upper curves), the characteristic magnetocon- 
ductance dip at zero magnetic field is seen. It is associ- 
ated with the effects of the coherent backscattering, or 
weak localization, which is present at B± = but is sup- 
pressed by a weak perpendicular magnetic field. Quanti- 
tatively similar <j(B±) dependences have been obtained 
on another sample, as shown in Fig. 2(b) (magnetocon- 
ductivity curves for this sample have been measured only 
at one direction of the magnetic field) . 

It is possible to fit the higher-n s magnetoconductance 
curves by the Hikami-Larkin formula [l£j 
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where r and are elastic and inelastic electron life- 
times, respectively, "P is the Digamma function, g v is 
the valley degeneracy, tb = H/4eB±D, and D is the 
diffusion coefficient; the numeric coefficient a depends 
on the ratio of intravalley to intervalley scattering rates 
and is expected to be between 0.5 and 1. An exam- 
ple of such fit is shown in Fig. 3 (a) (the top curve). 
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FIG. 1: Shubnikov-de Haas 

NH50.3064.9N at n s = 7.3 ■ 10 11 cm" 2 (a) and n s = 1.17 
10 11 cm" 2 (b). The vertical arrows in (b) show the positions 
of the SdH minima corresponding to Landau filling factors 
v = 14, 10, and 6 as labeled. Temperatures are indicated in 
the figures. Measuring current is 2 nA (a) and 200 pA (b). 
The experimentally measured amplitudes of the oscillations 
are shown in the insets by solid circles together with the cal- 
culated dependences (solid lines); the parameters used in cal- 
culations are indicated in the insets. The lower inset in (b) 
shows a close-up view of two SdH oscillations corresponding 
to v = 14 and 10 measured at three temperatures. 



The value of the phase-breaking time used, = 30 ps, 
is consistent with experimental data for (~ 20 to 
100 ps) obtained at similar temperatures on different 2D 
systems 0, 0, ^J, although it is more than an order 
of magnitude shorter than that expected theoretically, 
t « 27rft 2 cr/e 2 fc B Tln(crV2e 2 ) - 10~ 9 s [n|. 

As follows from Eq.^ the magnitude of the dip is equal 
to (e 2 ag v /irh) ln^/r). Using the above expression 
for the inelastic time, we obtain r^/r oc l/ln(er/i/2e 2 ). 
Therefore, the magnitude of the dip is expected to weakly 
(double-logarithmically) increase as the average conduc- 
tivity decreases, provided the variations in electron den- 
sity are small, as they are in our case. 

This is not what is observed in the experiment: as 
one approaches the transition, the magnitude of the dip 
sharply drops, and at the critical electron density (the 
lowest curves in Figs. 2 (a,b)), the dip is no longer seen 
on the scale of these figures. However, the shape of the 
magnetoconductivity does not change significantly with 
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FIG. 2: Longitudinal magnetoconductivity in a weak per- 
pendicular magnetic field at T — 42 niK for a range of elec- 
tron densities which are indicated near each curve in units of 
10 11 cm" 2 in sample NH50.3064.9N (a) and NH50.4064.7N 
(b). 



decreasing n s . This is illustrated by the middle curve 
in Fig. 3 (a) showing a(B±) multiplied by six (n s — 
0.90 • 10 11 cm" 2 ), which makes it quantitatively similar 
to the upper curve. This similarity demonstrates that 
the functional form of the <r(B±) dependence, described 
by the expression in brackets in Eq. ^ does not change 
noticeably as the density is reduced from 1.23 • 10 11 to 
0.90 • 10 11 cm" 2 ; instead, it is the vertical scale of the 
effect that rapidly decreases upon approaching to the 
MIT. At yet lower density, n s = 0.82 • 10 11 cm" 2 , the 
magnitude of the dip does not exceed 2% of that for 
n s = 1.23 • 10 11 cm" 2 (compare the upper and the lower 
curves in Fig. 3 (a)). Therefore, in the immediate vicinity 
of the metal-insulator transition, the coherent backscat- 
tering in the metallic phase appears to be strongly sup- 
pressed even in zero magnetic field . This is in agree- 
ment with the results obtained in ultra-high-quality p- 
GaAs/AlGaAs heterostructure [To|. 

In principle, suppression of the magnitude of the weak 
localization dip may be due to the inelastic scattering 
time being strongly reduced compared to that expected 
from the noninteracting theory. However, our attempts 
to approximate magnetoconductivity close to the tran- 
sition by the Hikami-Larkin formula with reduced 
have failed. An example of such an attempt is shown 
in Fig. 3 (b). The required magnitude of the dip can be 
achieved by choosing = 6 ps; however, the calculated 
dependence (the lower dashed curve) is much broader 
than the experimental one, and its curvature, d 2 u/dB 2 ^ 
is positive in the entire range of magnetic fields, contrary 
to the experiment. On the other hand, the experimen- 
tal curve in Fig. 3 (b) can be fitted reasonably well (the 
upper dashed curve) by choosing the same parameters 
as for the highest-density curve in Fig. 3 (a), but with 
the prefactor ten times smaller (a = 0.06). We therefore 
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FIG. 3: Magnetoconductivity in sample NH50.3064.9N at 
T = 42 mK and at different electron densities indicated near 
the curves in units of 10 11 cm -2 , (a) The curves are vertically 
shifted for clarity, and two lower curves are multiplied by 6 
(the middle one) and 50 (the lower one), which makes them 
quantitatively similar to the upper (unmodified) curve. The 
thick dashed line is the fit by Eq. Qwith a = 0.6 and = 
30 ps. (b) Unsuccessful attempt to fit the experimental data 
for n B — 0.85 ■ 10 11 cm" 2 by Eq. 0with the same prefactor 
a = 0.6 as above but reduced = 6 ps (the lower dashed 
curve, which is shifted for clarity). However, the data can 
be fitted reasonably well by choosing = 30 ps as in (a), 
but with the prefactor ten times smaller (a = 0.06; the upper 
dashed curve). 



believe that it is the value of the prefactor, rather than 
T0, which is responsible for the effect. 

As we have already mentioned, experimental data on 
Si MOSFETs reveal a sharp increase of the spin sus- 
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same region of densities where we observe the suppres- 
sion of the weak localization. Due to the strong renor- 
malization of the effective mass EES, the interaction 
parameter r s reaches at least 50 near n x (which in our 
samples practically coincides with the critical density for 
the metal-insulator transition). However, samples of a 
lower quality become strongly localized at densities no- 
ticeably higher than n x , where the effective mass re- 
mains close to its band value and r s does not exceed 10. 
Coherent backscattering in such samples is still present 
near the transition to a strongly localized state; there is 
no temperature-independent separatrix, and resistance 
of some of seemingly-metallic curves turns up at low 
temperatures |2l|. Similar behavior is observed in p- 
GaAs/AlGaAs heterostructures: "ordinary" weak local- 
ization corrections and low-temperature up-turns of the 
resistance have been reported in samples with relatively 
high transition density « 4.6 • 10 10 cm -2 |9j, while no up- 
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turn of the resistance and practically no weak localization 
corrections have been seen in higher-quality heterostruc- 
tures with almost an order of magnitude lower transition 
density « 6 • 10 9 cm- 2 10]. 

In summary, not too close to the transition, at 8 > 0.1, 
there are no principal deviations from the "conventional" 
magnetoconductance associated with the suppression of 
the coherent backscattering by a perpendicular magnetic 
field. This fact, although known from earlier measure- 
ments 0, 0, 0, |22, is not trivial, because the weak 
localization theory was developed for non- or weakly- 
interacting electron systems, while in our case, the energy 
of electron-electron interactions exceeds the Fermi energy 
by more than an order of magnitude. However, in the 



immediate vicinity of the MIT, where the role of interac- 
tions is especially strong, the magnitude of the weak lo- 
calization corrections sharply drops, and near the critical 
electron density, the corrections do not exceed 2% of their 
"noninteracting" value. This confirms similar results in 
an ultra-clean p-GaAs/AlGaAs heterostructure [Hj and 
thus establishes their universality. The absence of the 
localization corrections at the metal-insulator transition 
may account for the existence of a flat separatrix between 
metallic and insulating states observed in Refs.[T3, Hflj. 
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